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Abstract 

We present a study of the role of fermions in the decay of metastable 
states of a scalar field via bubble nucleation. We analyze both one and three- 
dimensional systems by using a gradient expansion for the calculation of the 
fermionic determinant. The results of the one-dimensional case are compared 
to the exact results of previous work. 



I. INTRODUCTION 

The study of the mechanism of decay of metastable systems via nucleation of thermally 
activated bubbles finds a wide range of applicability in Condensed Matter physics. In fact, 
the formalism developed in [1J and |2| fits naturally in the study of the dynamics of defects 
in one-dimensional conducting polymers []5]||. Applications to three-dimensional systems 
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are also numerous, from the behaviour of binary liquid mixtures to problems in Optics, 
such as atomic ionization in ultra-strong laser fields |I0| , and even Baryogenesis ||11|| . 



In a previous article 0, we studied a one-dimensional system of interacting fermions 
and bosons that started in a metastable vacuum and gradually decayed to the true one. 
We investigated the nucleation of bubbles of true vacuum inside the false one via thermal 
activation. We analyzed the stability of these bubbles and calculated the decay rate as a 
function of time, in the presence of a finite density of fermions, at finite temperature. The 
inclusion of fermions proved to have remarkable effects as it changed qualitative features of 
the physics of metastability p. 

The fact that we were working with a one-dimensional system allowed us to solve the 
problem exactly by means of inverse scattering methods. However, as we move to higher 
dimensions, this mathematical tool is no longer available and we are compelled to make 
approximations to evaluate the fermionic determinant that encodes the fermionic effects 
on the system. This may be accomplished by the use of a functional gradient expansion 
for the determinant. Throughout this work, we perform our calculations for both one and 
three-dimensional systems, so that we can compare the results of our approximations on 
the former with the exact results of [Q. In fact, our approximation is able to detect the 
phenomenon of "quantum stabilization" for both one and three-spatial dimensions. 

The paper is organized as follows: in section II, we present the model considered; in 
section III, we perform the gradient expansion for the fermionic determinant; in section IV, 
we obtain the form of the sphaleron solution that includes fermionic effects; in section V, we 
include the presence of the gap and two bound states of fermions; in section VI, we study the 
stability of bubble-like solutions; in section VII, we comment on the methods and present 
our conclusions. 

II. THE MODEL 

The model Lagrangian has the following form 
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£ = \{dM^<t>) - iv(<f>) - V{<k)] + Mh^ -n-gfl (2.1) 

where \i is the bare mass of the fermions, g is the coupling constant, ij) a {x) is the fermion 
field, a denotes fermion species, <f)(x) is a scalar field and <p2 is a local minimum of the 
potential 

V (<!>) = y(0-0o) 2 (V + 0o + y) + j4> (2.2) 

where 0o is a constant and j is an external current, responsible for the asymmetry of the 
potential even in the purely bosonic case. (We may find a physical realization of this form 
of potential in the description of conducting polymers see Fig. 1). 

We will be interested in the effects of fermions on the bosonic field. Thus, in order to con- 
struct an effective theory for bosons, we must integrate over fermions to obtain an effective 
action. However, integrating over fermions implies calculating the following determinant 

S F = -ln[det(i'fdf t - p - g<j>)] = -tr[ln{i^d^ -ft- g</>)] (2.3) 



In one spatial dimension, it is possible to perform an exact calculation of ( |2.3j ) by making 
use of inverse scattering methods F° r three spatial dimensions, we must resort to 

approximations. 

III. GRADIENT EXPANSION FOR THE FERMION DETERMINANT 

After the fermionic integration, we may rewrite the effective action as 



effl 



d"x ji(^)(«9^) - [V{<f>) - V(<h)]} - tr[ln(iYd, -ft- g(f>)) (3.1) 



where / d v x = Jq dt J dPx. The field configuration that extremizes (3.1) must satisfy the 
Euler-Lagrange equation 

dV 1 

□</> = -— -Sp<x\- — -\x> (3.2) 

d(f) -n-gp 

where Sp means trace over the spin degrees of freedom. 



We may calculate the Green function that appears in (3.2) by using a functional gradient 



expansion |T4|. This means that we will focus on the long distance (small momentum) 
properties of the theory. To do so, we use the identity fll"2"| , |T3f (see Appendix A) 
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Q(x, x) = Sp < x 
= Sp 

where 



\x > 



i^dy. + M(x) 
d v r> 1 00 ( 1 8 
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(27r)" 7 ^ + M(a;)^ v 



i dp j^Pfj, + M(x) / 



(3.3) 



AM(~,i) = d,M(x)-^- + l -d,d u M{x){-) 2 °~ 



i dp' 



i dp. 



i dp il dp v 



(3.4) 



Keeping terms up to second order in the derivatives and explicitly performing the inte- 
grals, we obtain 



G(x, x) = a u M u ~ L + %M V -\/ 
where a v and (3 V are functions of the number of space-time dimensions: 



(3.5) 



a,, 



2 u+l 



■r(i - 1//2) 



(3.6) 
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-i//2 



2^- 1 r(z//2) 



w>. r(4 _ „ /2) _ ni±iV!) r( 3 _ „ /2) + £<i±^ r(2 _ „ /2) 



4! 



2!2! 



4! 



(3.7) 



For v — 2 and z/ = 4, these functions are divergent. Nevertheless, we may absorb 
these divergences by a suitable redefinition of the free parameters of the model. After a 
redefinition of g, \x and j in (2.2) and remembering that, in our case, M(x) = —(// + g<p), 
we may explicitly calculate the equation of motion (3.2). 

For v = 2, equation (3.2) reads (up to order g 2 ) 



<9V 



(/i + #) 2 



(A* + 9<P) 2 + £ 



(3.8) 



and, for v = 4, 



The part of the effective action, Sp, associated with the fermionic determinant has the 
following form |TJJ 



S F = -g J d u x J [D<j>\ Q{x,x) (3.10) 
with the matrix element given by (3.5). Thus, for v = 2, 

S^-J-lfx-r^, (3.11) 

We note that, in the case of u = A, the only effect of the fermions is (to this order of 
approximation) to renormalize the free parameters of the theory. We will see, in section 
V, that this situation is completely changed if we include a gap and bound states in the 
fermionic spectrum (see Fig. 4). 

IV. THE SPHALERON SOLUTION 

Based on the results of our previous work, we shall look for sphaleron-like solutions to 
the equations of motion obtained in the last section. A sphaleron-like solution corresponds 
to a field configuration that starts at the false vacuum 02, almost reaches the true vacuum 
0i, and returns to 02, i.e., it looks like a droplet [|l],0,|6],|9|] . 

A. The v = 2 case 

In this case, the equation of motion is (3.8). Defining ip = + /VflS using the thin- 
wall approximation [|[ (i.e., considering a droplet whose radius is much larger than its wall 
thickness) and imposing the boundary conditions 

<-Ps P h(x -> ±oo) -> ip 2 (4.1) 

d(fsph 



dx 



(x -> ±oo) -> (4.2) 



we obtain 



Vsph = ¥2- <t>p[tanh(£ + f ) - tanh(£ - f )] 



(4.3) 



where 



1/2 



(4.4) 



i = g(p P (x - x c . m ) 



(4.5) 



/ 



6 = -cosh 1 



(4.6) 



\ 1 

) 

and where we have assumed /i/g » 1 (neglecting corrections of order O ((g / /i) 2 )) , consistent 
with the conditions of validity of the gradient expansion, in order to obtain a closed form 
for the function that represents the sphaleron. 

The parameter x c _ m , reflects the translational invariance of the equation of motion and £0 
is related to what is usually called the radius of the sphaleron (see Fig. 2), being extremely 
important in the analysis of stability (see section VI). 



B. The v = 4 case 

In this case, the equation of motion is (3.9). Thus, assuming a solution with radial 
symmetry, we have, in the thin-wall approximation, 



fsphii) =¥2- <j>p[tanh(£ + |o) - tanh(i - £ )], I > 
where we have imposed the following conditions 



(4.7) 



(4.8) 



(4.9) 
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and 



Vsphij- -> oo) -> y? 2 (4-11) 

_ oo) _ (4.12) 

Therefore, we have a true vacuum bubble, of radius £ > centered at the origin (see Fig. 
3). The situation is the same as we would have encountered in a purely bosonic system [|TJ. 

V. INCLUSION OF A GAP AND BOUND STATES 

In the previous section, we have shown that a bubble-like configuration still satisfies the 
modified equations of motion (In fact, in the v = 4 case, we have no changes at all). These 
bubbles play the role of a background field in the determination of the fermionic spectrum. 
In fact, in some cases, instead of a simple continuum, we may have the presence of a gap and 
pairs of bound states (we will consider, for simplicity, just one pair). Therefore, if we 
have a richer spectrum (see Fig. 4), the calculation of the fermionic determinant, performed 
in section III, should be reviewed to deal explicitly with the bound states. The calculation 
that led us to (3.11) is the result of a gradient expansion approximation for the total trace 
(2.3), where we sum over all continuum and bound states. However, as we are interested in 
an effect that depends crucially on the relative occupation of the bound states H , we should 
consider a finite density of fermions, i.e., a partial trace. Thus, instead of summing over 
all fermionic momenta (the restriction of small momenta is over the bosonic fields only), we 
shall sum them only up to the (occupied) top bound state (see Fig. 4(c)). 

Therefore, the complete effective action has the following form 

SeffW = J d?x - [V(<f>) - V(<h)]} - 

- 9 -jd v xj[D4>] G(x,x)- 



- n +9 J Q dh J dx dt 2 J dy J [D<j>\ ^* B (x,ti, [4> sp h\) x 

x g(x, y)ip B (y, *2, [(psph]) (5.1) 

where the first term is the bosonic contribution, the second term represents the "Dirac sea" 
(including antiparticle bound states) and the last one is due to the particle bound states. 
TpB is the wavefunction of a bound state and n + is its occupation number ("doping"). 

The only term that remains to be calculated is the last one. Assuming that the occupation 
of the bound states will not affect in an appreciable way the form of the bubble (in the one- 
dimensional case, this is an exact result 00), we may rewrite this term as (see Appendix 
B) 

Sbound = n + gT J dxp{x)(f) sph (x) (5.2) 

where p(x) = i[) b (x)i/j b(%) is the normalized density of probability distribution of the bound 
charge. However, it is already known that the charge associated with a bubble tends 
to concentrate on its surface in a gaussian-like way. For our purposes, we will assume that 
a delta-like distribution will be a reasonable approximation. 



A. The v = 2 case 



In this case, we may write the density p as 
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(5.3) 



The energy of the bubble, E = —S e ff/T, as a function of the radius £o, has the following 
form: 



£(£o 



"+°° d£ J 1 
-oo g(j) P 1 2 

— n + tanh(2^o) 



1 + 



1 9 



12vr \h 



+ iv(<p sph ) - V(<h)] 



(5.4) 



where we have incorporated the contribution of the Dirac sea in the first term, by making use 
of the equation of motion. In fact, the contribution of the continuum is of order 0((g/p) 2 ) 
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and may be neglected to this order of approximation. The only relevant contribution of 
fermions comes from the bound states. 



B. The v = 4 case 



In this case, we may write the density p as 



\2 



P = ^ffH(£ - £o) (5-5) 

The energy of the bubble as a function of the radius £o has the following form: 

2 



B(&)M f i^'nr) 

— n + tanh{2£ > o) (5.6) 



VI. STABILITY OF THE BUBBLES 

To analyze the stability of the bubble-like solutions obtained in section IV, we shall 
study the behaviour of the the energy of the bubble as a function of the bubble radius £ > 
now considered as a dynamical variable, s. The results, for the cases v = 2 and v = 4 
obtained in last section, are plotted in Fig. 5 and Fig. 6, respectively. Our results for 
the v — 2 case should be compared with the exact results, previously obtained by using 
inverse scattering methods (see Figs. 5(b) and 5(c)). In this way, we may control our 
approximations. In fact, the observation of Fig. 5 shows that our approximation preserves 
the "quantum stabilization" brought about by fermions. Nevertheless, we find a quantitative 
difference between the approximate and the exact results due, mainly, to our naive delta-like 
approximation for the fermionic density. In order to improve on this approximation, one 
may use a gaussian-like pattern or even the exact fermionic density in the presence of the 
sphaleron. However, in doing so, one can no longer use the simple analytic form of the 
previous section. 
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VII. CONCLUSIONS AND COMMENTS 



The aim of this work was to understand the role played by fermions in the decay of 
metastable states of a scalar field via bubble nucleation. Our basic approximation was a 
gradient expansion for the fermionic determinant that proved to work well in a comparison 
with exact results previously obtained for one-dimensional system. From our results, it is 
clear that the effects of fermions in an arbitrary number of space dimensions are due, al- 
most exclusively, to the relative occupation of bound states. The states of the continuum 
contribute only to order 0((g/ [i) 2 ), and may be neglected. (This is in agreement with re- 
sults obtained for an analogous problem in QCD ||15||.). For v = 2 and v = 4, the new 
charged bubbles have the same functional form as the purely bosonic ones, except for a 
reparametrization, and have their stability drastically modified by the fermions: besides un- 
stable bubbles, we find metastable bubbles as a result of a "quantum stabilization" brought 
about by the fermions. For v = 2, this result was obtained exactly in Ref. 0. This strik- 
ing feature may, in principle, be measured in some realistic systems. The two-dimensional 
case finds a natural application in the physics of linearly conducting polymers 0-0]. The 
results for the four- dimensional case may have important consequences for the physics of 
Baryogenesis and in some problems in Optics [P!D] , PH| 

From the results for the energy of the bubbles as a function of their radii, it is possible 
to calculate decay rates for the metastable states of the scalar field as explicit functions of 
time. This may be implemented by using the formalism presented in Q and and allows 
for a non-equilibrium description of the decay process for both v = 2 and v = 4 cases. 
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APPENDIX A: GRADIENT EXPANSION OF THE FERMION DENSITY 



The complete derivation of identity (3J3) can be found in ref. [12]. For completeness, we 
present here the main steps of its demonstration. 

The density of fermions at a given point xo, in a background field M(x), has the following 
form 

p(xo) = g(xo,xo) = Sp<x \——^———\x > (Al) 

27 M Cyj + M [x) 

where |x > is a position eingenstate with eigenvalue x , and Sp is the trace over spin degrees 
of freedom. 

In momentum representation, this expression takes the form 

P[Xo) * P J (2ny e YP» + M{x) [ ] 

In order to transfer the Xq dependence to M(x), we make a unitary transformation that 
eliminates the exponentials from the last expression: 

P(x " ) - Sp l^ 1 ^ + M(i + X0} (A3) 

Note that Xq in M(x + Xq) is a c-number, not an operator. 
We can now expand M(x + x ) around x : 

M(x + x ) = M(x ) + AM(x ,x) (A4) 

where 



AM(x , x) = d^M(x Q )x^ + -d li d v M{x Q )x^x v + ... (A5) 



Using (|A4|) in (|A3|) and considering that, in momentum space, x^ = we have 

r d v p 1 
P(X0) = SP J (2 7 r)^^ + M(xo) + AM(x ,i| ; ) {M) 

We can factor out —n — hnr, — T obtaining: 
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p(x ) = Sp 



d v p 



(27t)-7^ + M(x ) 



1 d 

1 + AM(x ,-— ) x 



(A7) 



The idea is to expand the last parenthesis in powers of AM/ (7 M p M + M) . This means 
that we consider smooth backgrounds M(xq), so that \d jjl M(xo)\ < \M(xq)\. In momentum 
space, this condition implies \k^M(k)\ < \M(k)\, where k < 1 is the momentum transferred 
from the fermions to the boson background. With this expansion we finally obtain: 

p(x ) = Sp [ -4 tvt-^ Y ("l) n AM(x , --%-) x — - — (A8) 

APPENDIX B: ENERGY OF BOUND STATES 



In this appendix we derive equation (|5.2|) . 
From (5.1) we have 

Sbound = -n+g J dh J dx J dt 2 J dy J [D(p] x 

x ^b(z, h, [<P v h])G(x, y)Tp B (y, t 2 , [<psph}) (Bl) 

where ipB is the bound-state wave function in the presence of a sphaleron, and can be written 
in the following form 

M^MM) = M^ iEBt ( B2 ) 

where Eb is the energy of the bound state. 



The Green function in ( |B1| ) takes the form 

g(x, y)=J2 ^ n (y)Mx )e~^-^ (B3) 

n 

where E n < ep and 

(iffy + fj, + (f))M%) = EnM*) (B4) 

The index n may be discrete or continuous, depending on the part of the spectrum we are 
considering. The sum must be done up to the Fermi energy e_p. 
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Using (p3|) in (pl|), we obtain 



Sbound = -n+g J dx J dy J dt x J dt 2 J [D<j>\ 



x 

n<tp 

x ^ B (x)M^i(y)My) e-^-^-^ (B5) 

The difficulty in evaluating this expression resides in the ^-dependence of ip n . However, 
in the semiclassical context we consider fields "near" the sphaleron configuration. We also 
assume that the occupation of the bound states will not affect essentially the form of the 
bubble (this is an exact result in one spatial dimension 0). Considering static fluctuations 
around the sphaleron, we can interchange the order of evaluation of the time and the 
integrals and, using the formula 



/ dt x J dt 2 e-^-^K^) = TS ^ En _ Eb) 
we obtain 

Sbound = -n+gT J dx J dy J [D<f>] ^\{x)^ B {x)^ B {y)^ B {y) (B6) 

provided that €f > Eb- 

Note that the S(E n — Eb) fixes the ^-dependence of if) n , so that we can integrate in <fi 
and then evaluate the expression in the sphaleron configuration. This procedure leads to 

Sbound = -n+gT dx dy (j) sph {x^^ B {x^^ B {^ B {^ B {y) (B7) 



In this form, the x and y integrations decouple and, using orthonormality of the wave 
functions, we finally arrive at 

Sbound = -n+gT J dx (j) sph {x)ip B (x)il} B (x) 

= -n + gT J dx p(x)(j) sph (x) (B8) 



which is precisely equation (|5.2| ). 
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Figure Captions: 

Figure 1: Form of the potential V{4>). 
Figure 2: The "sphaleron" solution for v — 2. 
Figure 3: The "sphaleron" solution for v — 4. 

Figure 4: (a) Continuum (gapless) fermionic spectrum, (b) Fermionic spectrum in the 
presence of a constant background, (c) Fermionic spectrum in the presence of a bubble-like 
background. Q is a cutoff (in polymer applications, it corresponds to the bandwidth) 

Figure 5: Energy of the bubble as a function of the bubble radius for v — 2. (a) purely 
bosonic. (b) zoom of the interesting region for the exact solution (including bound states), 
(c) zoom of the interesting region for the approximate solution (including bound states). 

Figure 6: Energy of the bubble as a function of the bubble radius for v — 4. (a) purely 
bosonic. (b) zoom of the interesting region for the approximate solution (including bound 

states). 
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Figure 5(b)-exact 




Figure 5(c)-approximate 
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